Real-time diagrammatic approach to transport through interacting quantum dots 

with normal and superconducting leads 
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We present a real-time diagrammatic theory for transport through interacting quantum dots 
tunnel coupled to normal and superconducting leads. Our formulation describes both the equilibrium 
and non-equilibrium superconducting proximity effect in a quantum dot. We study a three-terminal 
transistor geometry, consisting of a single-level quantum dot tunnel coupled to two phase-biased 
superconducting leads and one voltage-biased normal lead. We compute both the Josephson current 
between the two superconductors and the Andreev current in the normal lead, and analyze their 
switching on and off as well as transitions between 0- and 7r-states as a function of gate and bias 
voltage. For the limit of large superconducting gaps in the leads, we describe the formation of 
Andreev bound states within an exact resummation of all orders in the tunnel coupling to the 
superconducting leads, and discuss their signature in the non-equilibrium Josephson- and Andreev- 
current and the quantum-dot charge. 

PACS numbers: 74.45.+c,73.23.Hk,73.63.Kv,73.21.La 



I. INTRODUCTION 

Continous advancements in nanofabrication have made 
it possible to attach superconducting leads to quan- 
tum dots. The supercurrent through a quantum dot 
has been measured through dots realized in carbon 
nanotubesi and in In As nanowires* 2 - Recently, transport 
measurements on a single self-assembled InAs quantum 
dot coupled to Al superconducting electrodes have been 
reported^ 

From a theoretical point of view, quantum dots cou- 
pled to superconducting leads are of great interest, since 
a very rich physics is expected from the combination 
of superconducting correlations, electron-electron inter- 
action and non-equilibrium in the dot. Subgap trans- 
port through a normal-dot-superconductor system is sus- 
tained by Andreev reflection 4iai&Z&£ The Josephson 
coupling between two superconductors through a quan- 
tum dot has been addressed in the limit of a non- 
interacting quantum dot in Ref. 110 . In the opposite limit 
of a large charging energy, the electrons forming a Cooper 
pair tunnel one by one via virtual dot states ; 11 ' 12 : 13 which 
establishes a Josephson current carried by higher-order 
tunneling processes. Other aspects of the problem, such 
as the Kondo regim o 11 : 14 : 15 : 16 : 17 : 18 : 19 : 20 or multiple An- 
dreev reflectio n 21 i 22 have also been addressed. The de- 
pendence of the charge in the quantum dot on the gate 
voltage and on the superconducting phase difference has 
been investigated in Ref. Moreover, numerical ap- 
proaches based on the non-crossing approximation^ the 
numerical renormalization group^ and Monte Carlo^S 
have been employed to study trans por t through this type 
of systems. The authors of Ref. [27] compare different 
approximation schemes, such as mean field and second- 
order perturbation in the Coulomb interaction. The 
proximity effect in double-dot systems has been also in- 



vestigated in different regimes . 17 i 18 i 28 In Ref. the non- 
equilibrium Josephson and Andreev currents through a 
dot coupled to one normal and two superconducting leads 
have been studied in the weak-proximity limit, consider- 
ing only first-order processes in the tunnel coupling with 
the superconductors. In this regime, finite Josephson and 
Andreev currents can flow only if the dot is driven out of 
equilibrium. The idea of using non-equilibrium to control 
the behavior of a Josephson junction has been proposed 30 
and experimentally tested 3 ^ some years ago. 



In the present work, we develop a real-time trans- 
port theory for an interacting quantum dot connected 
to both superconducting and normal leads. The theory 
can be conveniently formulated by means of a diagram- 
matic language and it is suitable for dealing with su- 
perconducting correlations, strong Coulomb interaction 
and non-equilibrium due to arbitrary bias voltages on 
the same footing. We demonstrate the use of our formal- 
ism for two examples. First, we study the equilibrium 
Josephson current between two superconductors due to 
cotunneling through the quantum dot and analyze the 
formation of a 7r-state for increasing on-site Coulomb re- 
pulsion on the dot. Second, we consider a transistor ge- 
ometry with one normal and two superconducting leads 
with large superconducting gaps. We calculate the non- 
equilibrium Josephson and Andreev current to all orders 
in the coupling strength with the superconductors, where 
the quantum dot is driven out of equilibrium by applying 
a bias voltage to the normal lead. This geometry is suit- 
able for performing a spectroscopy of the Andreev bound 
states in the interacting quantum dot. 
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II. MODEL AND FORMALISM 
A. Hamiltonian 

We consider a single-level quantum dot tunnel coupled 
to both normal and superconducting leads. The total 
Hamiltonian of the system is given by 



H = Hr 



■H t , 



(1) 



N is the total number of electrons, which equals the 
number of Bogoliubov quasiparticles plus twice the 
number of Cooper pairs. In the case that 77 refers to a 
normal lead, the order parameter vanishes, A,, = 0. 

The coupling between the dot and the leads is taken 
into account by the tunneling Hamiltonians 



H 



tunii,?? 



V, 



fe<7 



H.C. 



(7) 



The different (superconducting or normal) leads are la- 
beled by the index 77. The quantum dot is described by 
the Hamiltonian of the single-level Anderson model, 



Hjy = y~] ed) a d a + Un^ni 



where d a is the annihilation (creation) operator for 
an electron in the dot, e denotes the energy of the single- 
particle level, n a — d\d a is the number operator for spin 
a =t, 4, and U is the energy cost for double occupation. 
The leads' electrons are described by the annihilation and 
creation operators c^kcr and cL ka , respectively. In addi- 
tion to the kinetic-energy term ^2 ktJ e kcl lk(7 c rjka m the 
Hamiltonian, there may be a BCS pair-interaction part 

~9vJ2k,k' c l,k-\ c \-ki c v~k' ic v k'i to account for supercon- 
ductivity. On the one hand, we want to treat the inter- 
action on a mean-field level. On the other hand, we want 
to keep track of the total number of electrons, which will 
be important for situations with finite bias voltage be- 
tween different superconductors. This can be achieved 
by representing the lead electrons in terms of Bogoli- 
ubov quasiparticle operators 7^? CT and Cooper-pair an- 
nihilation (creation) operators^ via the Bogoliubov 
transform 



7r;fcT 



with coefficients 



V rjk :y r} u rjk 



\ 




V(efc-^) 2 + |A,| : 



v k = e 



8*7, 

\ 



2 I V^-^ + IA,,! 2 



(3) 



(4) 



(5) 



where fj, v is the electrochemical potential of lead 77 
and $ r) the phase of the order parameter A„ = 
9riJ2k(^l c v-kl c nk'\)- As a result, the mean-field Hamil- 
tonian for lead 77 reads 



(6) 



where for the sake of simplicity the tunnel matrix ele- 
ments Vfj are considered to be spin and wavevector in- 
dependent. The tunnel- coupling strengths are defined as 



(2) r r) = 2ir\V v \~ J2k < ^( a -'~ e fe)j which we assume to be energy 



plus an irrelevant constant. Here, E vk = of the re 
\J{ck — Uri) 2 + IA,,! 2 is the quasiparticle energy, and P$ = P|. 



independent. 



B. Diagrammatic Real-Time Technique 

The main idea of the diagrammatic real-time technique 
is to integrate out all the (noninteracting) fermionic de- 
grees of freedom in the leads to arrive at an effective de- 
scription for the reduced system, that is characterized by 
the state of the quantum dot and the number of Cooper 
pairs in the superconducting leads. The Hilbert space 
of the single-level quantum dot is four dimensional: the 
dot can be empty, singly occupied with a spin-up or spin- 
down electron, or doubly occupied. These are denoted by 
\x) € {|0>, I T),| i),|D) = djd\\0)}, and have energies E , 
= Ei, and Ed, respectively. The condensates in the 
superconducting leads are characterized by the number 
of Cooper pairs |n), relative to some arbitrarily chosen 
reference, where n is the vector of Cooper-pair numbers 
n v for each superconducting lead 77. The energy con- 
tribution from the Cooper-pair condensates is given by 
E n = J^ri^v^v ^ a ^ superconducting leads are kept 
at the same chemical potential then this energy contribu- 
tion simply provides a trivial additive constant. For finite 
bias voltage between at least two superconducting leads, 
however, the total energy depends on how the Cooper 
pairs are distributed among the superconducting leads. 

We start with the full density matrix of the total sys- 
tem, including the quantum dot, the fermionic degrees 
of freedom of the leads, and the Cooper-pair conden- 
sates. Since the fermionic degrees of freedom in the 
leads act as reservoirs, we can trace them out to obtain 
the reduced density matrix p ro( j with matrix elements 

Here, the label £ = (x, n) with energy E^ = E x + E n 
includes both the quantum-dot state x and the number of 
Cooper pairs, n, in the leads. For the diagonal elements 
of the reduced density matrix we also use the notation 
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1. Kinetic equation and current formula 

The dynamics of the reduced density matrix is gov- 
erned by the kinetic or generalized master equation, 



dt ^ w 



E 



^2 



(8) 



where the kernels Wv 1 !, 1 ft, t'l describe transitions due to 
tunneling. The current in lead rj can be written as 



J v {t) 



- e E 



dt'wgt n {t,t')Plt{t'), (9) 



is the sum of all kernels that describe transitions in which 
in total s electrons are removed from lead rj. 

Both the generalized master equation and the expres- 
sion for the current can be further simplified when all 
voltages and coupling strengths are kept time indepen- 
dent. The kernels do, then, only depend on the time 
difference t — t', and to determine the DC component of 
the current and all density matrix elements we only need 
the time integrals of the kernels, which we refer to as gen- 
eralized rates wH'/fi = f* dt'Wl^Ht - t') and gen- 

cralized current rates wff^ /h = dt'W^ v {t - t'). 

The indices £ contain more information than needed 
for our purpose. This is related to the fact that only 
the change and not the absolute value of the number of 
Cooper pairs in each superconducting lead matters, i.e. 

the value of the generalized rate If, t , does not change 
when we perform the simultaneous shift ni — > ni + m, 
112 — > n2 + m, n'i — > n'i + m, and n'2 — > n'2 + m for a 
given vector of additional Cooper-pair numbers m. After 
defining 



P£(m) ee £P ( 



(X2,n) 



W^Mm, m') 



^(X2,n)(x 2 



- .-n)(xi,m'+n') 
(X2,n)(x 2 ,n') 



(10) 

(11) 



and similarly for the generalized current rates, we obtain 
for the stationary current in lead rj 



E 

xxix 2 n ' 



W 



XXiV 
XX 2 



(0,n')P x /(n'), 



(12) 



where the matrix elements P v , 1 (n / ) are determined from 
i(E Xl -E Xa + E n )Pg(n) = J2 ^x X 2x1( n ' n Xl( n ') 



xix 2 n ' 



(13) 



together with the normalization condition Ey-^xW = 
1. Note that, in Eqs. (TT2"|) and (TT^J), due to conservation 
of the total number of electrons, only those Cooper-pair- 
number vectors n appear for which ^2 n v equals twice 
the number of dot electrons in state \i minus that in 
state xi (and the same holds true for n', X2: and Xi)- The 
generalized master equations for Pg* (n) with all other 
vectors n, not satisfying the condition stated above, de- 
couple and are, therefore, irrelevant. For illustration, let 
us consider the matrix element P(P(n) in a system with 
two superconducting leads; for example, with n = (—1, 0) 
or n = (—2, 1) it contributes, while with n = (1,0) it is 
irrelevant. 

In the special case that all superconducting leads are 
at the same chemical potential jis, the situation simpli- 
fies further. Due to the fact that in Eq. (IIB"]) the en- 
ergy contribution P n is the same for all n that are com- 
patible with xi and \2, the generalized master equation 
Eq. ([TB")) remains unchanged under the shift n — > n + m 
and n' — > n' + m with Er/ TO »? = 0- As a consequence, 

the generalized rates W XlXl and W xx } n as well as the 

° X2X2 xx 2 

solution for P^ 1 become independent of the Cooper-pair 
numbers, i.e., we can simply drop the arguments 0, n 
and n'. It is this limit that we are going to analyze in 
the results section of this paper. 



2. Time evolution of the reduced density matrix 



We generalize the real-time diagrammatic approach to 
transport through interacting quantum dots of Ref. [U 
to the case of superconducting leads. Our goal is to give 
a diagrammatic prescription to compute the generalized 

rates W XlX / (n, n') and W xx /"(0,n'). For this, we an- 

X2X2 v ' ; xx 2 v ' ; 

alyze the time evolution of the reduced density matrix 
that we obtain by integrating out the fermionic degrees 
of freedom in the leads. 

We assume at some initial time to (with to — > —00) 
the total system to be in a product state of the leads' 
fermionic degrees of freedom (taken at equilibrium) and 
the degrees of freedom of the reduced system. The 
time evolution of the reduced density matrix from time 
to to time t can, then, be described by Pj^ 1 (t) = 

E£ ic j ngg (t, t )Pf (to)- The propagator (t, t ) 

can be computed by means of a perturbation expansion 
in the tunneling Hamiltonian Ptunn = Ptunn,?;- How 
this is done has been explained elsewhere^ and here 
we will limit ourselves to sketch briefly the derivation, 
thereby pointing out the new ingredients due to super- 
conductivity. The propagator (starting from and ending 
at a product state of the leads' fermions and the reduced 
system) is written in interaction representation with re- 
spect to Ptunn as 



4 



n|§ (t, = Tr lcads { (& (f )|T K [|6 (f )) (6 (*)|e 



4/ dt"/f tunn (t")i 



(14) 



r 



being K t /_ 1 . i the Keldysh contour going from t' to £ and 
then backwards to t', Tk the time-ordering operator on 
the Keldysh contour, H tunn (t)i the tunnel Hamiltonian 
in interaction representation, and Tri ea( j s the trace over 
the fermionic part of the lead degrees of freedom. 

Next, we expand the exponential function in a power 
series of -fftunn- Finally, we perform the trace over the 
fermionic lead degrees of freedom by means of Wick's 
theorem. This is possible because the Hamiltonians of 
the leads are quadratic in the lead fermionic operators 
(this applies also to the superconductors in the mean- 
field description adopted here). 

For normal leads, only contractions between electron 
creation and annihilation operators are non zero. They 
are graphically depicted as tunneling lines (normal lines) 
with an arrow going from the vertex oL kg .d a to the vertex 

dl.c v ka- We define the direction of the line such that an 
electron is removed from the dot at the vertex where the 
line starts and added to the dot at the vertex where the 
line ends. 

For the superconducting leads two different types of 
lines appear. There are again normal lines, connect- 
with the dJ.Cjjfco- . In addition, 



mg a vertex c^ ka d. 



there are anomalous lines, connecting either c^ fe(T d 



C l-k-a d - 



dt 



'■rjka 



with d 



~"q—k- 



with 
Due to the con- 



vention for the arrow direction introduced above, the 
anomalous lines carry two arrows that point towards each 
other (outgoing anomalous line) if two annihilation op- 
erators of dot electrons are involved, and away from each 
other (incoming anomalous line) for two creation opera- 
tors of dot electrons. To evaluate the contractions, it is 
convenient to perform the Bogoliubov transform for the 
lead electron operators. Since the Bogoliubov transform 
involves the operators , the number of Cooper pairs 
in lead n may be changed at the tunnel vertices. For each 
normal and each anomalous superconducting line, there 
are two possibilities. The vertices being connected by 
normal lines either involve no operator or one S n and 
one SJ.. For outgoing (incoming) anomalous lines, either 
one of the two vertices carries the operator Si (S v ). The 
normal lines describe quasiparticlc tunneling. Whereas, 
the anomalous lines describe Andreev tunneling: for an 
incoming (outgoing) anomalous line a Cooper pair breaks 
(forms) in the lead and its constituents enter (leave) the 
dot at the two vertices. 

Now we can describe the different contributions to the 
propagator by means of a graphical representation on the 
Keldysh contour. An example is shown in Fig. [TJ where 
both normal and anomalous lines are present. 

The propagator obeys the Dyson equation 



J 



n||(t,t') = n||W(t,tO% el ^ 



2«2 



E 

tig 



dt" / dt"'U 



66 



(15) 



where n! 1 ! 1 ^ (t, t') is the free propagator and we identify only keep the information of Cooper-pair-number diffcr- 



"«2 6 

W tltl (*' *') with thc irreducible 

part of the propagator, 

i.e. with the sum of irreducible diagrams going from t' 

to t contributing to the propagator U^}(t,t'), where a 

diagram is irreducible if any vertical line through the 
diagrams cuts at least one tunneling line. The order in 
r of an irreducible diagram is given by the number of 
tunneling lines present in the diagram. Examples of first- 
and second-order diagrams are shown in Fig. [1] 

From now on we concentrate on stationary situ- 
ations and, hence, we will consider the generalized 
rates. These are given by the Laplace transform of 



W^(t - *') computed 



at z = 0+ 



IT"' 



6Ci _ 



r* dt'e-^-^wfikt-f) 

■' ^ ?2?2 v ' 



«2?2 

Furthermore, we 



ences, i.e., we formulate the rules for W ,(m, m') = 

X2X2 v ' 

En n' ^ (xi '™ + ? )( ^' m ' +n ' ) . The last step does not only 

^ n > n (X2,n)(x2,n) 1 J 

reduce the number of matrix elements to be considered. 
Another virtue is the possibility to combine different 
contributions. As mentioned above, for a given super- 
conducting line there are always two possibilities to as- 
sign operators to the two vertices. Depending on 
the topology of the reduced diagram, the correspond- 
ing terms may contribute to different generalized rates 

wl^/ 1 (that differ from each other by the number of 
Cooper pairs) but they always contribute to the same 
W X1X } (m, m'). It turns out that, since the tunneling 
strengths T n are independent of energy, the analytic ex- 



z=0+ 



■5 
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AEi AE 2 AE 3 AE 4 AE b 

FIG. 1: Graphic representation of a contribution to the el- 
ement iljij of the propagator for the reduced density matrix, 
for the exemplary case of one superconductor with chemical 
potential /is. The upper and the lower line of the Keldysh 
contour represents the forward and the backward propaga- 
tion, respectively. From left to the right we can identify a 
first-order diagram with a normal line, a first-order diagram 
with an anomalous line, a second-order diagram with a nor- 
mal and an anomalous lines. Below the diagrammatic repre- 
sentation of the propagator, we show our arbitrary choice for 
the directions of the anomalous lines and the corresponding 
energy differences to be used in rule 2, which are given by 
A£i = wi - e; AE 2 = -0J2 - e + 2/* s ; AE 3 = -w 3 - e + 2/x s ; 
AE4 — 0J4 — 0J3 — 2e — U + 2/is ; AE 5 = CJ4 — e. 



pressions of the two contributions always combine nicely, 
which leads to a rather compact formulation of the dia- 
grammatic rules presented below. 

As defined above, the arrows of a tunneling line indi- 
cate whether an electron enters or leaves the dot at a 
given tunnel vertex. Furthermore, we want to define an 
overall direction of each tunneling in order to define the 
sign of the energy carried by the Bogoliubov quasiparti- 
cles. For normal lines, we will always choose the direction 
set by the single arrow. For anomalous lines (that carry 
two opposite arrows), we pick the direction arbitrarily, 
and assign the creation or annihilation of a Cooper pair 
to the vertex at which the line direction is opposite to 
the arrow. 

In order to construct a systematic perturbation expan- 
sion in the tunnel coupling, both the generalized rates 
and the probabilities are expanded in orders of T, i.e 
(m, m') = (1) (m, m') + (2) (m, m') + 

G(r 3 ) and P£(m) = P x x 2 l(0) (m) + + 0(T*), 

where the superscript indicates the order in T. 

3. Diagrammatic rules 

The rules for evaluating the generalized rates 

W XlX 'Hm,m') are: 

X2X2 v ' ; 

1. Draw all topologically different diagrams with fixed 
ordering of the vertices in the real axis. The ver- 
tices are connected in pairs by tunneling lines car- 
rying energy Wj. The tunneling lines can be normal 



or anomalous. For each anomalous line choose the 
direction (forward or backward with respect to the 
Keldysh contour) arbitrarily. 

2. For each vertical cut between two vertices assign 
a factor 1/(AE + it]) with r] = + , where AE 
is the difference between the left-going and the 
right-going energies, including the energy of the dot 
states, E x , the tunneling lines, u>i, and the energy 
difference of Cooper-pair condensates, E n . The lat- 
ter is increased (decreased) at each vertex of an 
outgoing (incoming) anomalous line at which the 
arrow is opposite to the arbitrarily chosen line di- 
rection. 

3. For each tunneling line assign a factor 
^7 1 D v (uJ l )f^(uJ l ), where /+(w t ) = f v (ui) = 
[l+cxp(a) i -/i r ,)/(fcBT)] _1 and/~(a>i) = l-f v (u>i), 

and A,M - Tsdgfci^ l" " ^ - l A "»- 
The upper (lower) sign applies for lines going 
backward (forward) with respect to the Keldysh 
contour. For anomalous lines multiply an ad- 
ditional factor ±sign(cjj)^ 2 |. Moreover, assign 

a factor e -1 *" for an outgoing and e 1 *" for an 
incoming anomalous line. [For normal leads, only 
normal lines with D v (u>i) = 1 appear.] 

4. Assign an overall prefactor —i. 
Furthermore, assign a factor —1 for each 

a) vertex on the lower propagator; 

b) crossing of tunneling lines; 

c) vertex that connects the doubly occupied dot 
state, |D) = rf|<i{|0), to spin up, | f); 

d) outgoing (incoming) anomalous tunneling line in 
which the earlier (later) tunnel vertex with respect 
to the Keldysh contour involves a spin up dot elec- 
tron. 

[The factors in c) and d) arise due to Fermi statis- 
tics from the order of the dot and lead operators, 
respectively.] 

5. For each diagram, integrate over all energies a;,. 
Sum over all diagrams. 

The generalized current rates W xx , lV (0, m') are evalu- 
ated in the following way: 

6. Multiply the value of the corresponding generalized 
rate W**, 1 (0, m') with a factor given by adding up 
the following numbers for each tunneling line that 
is associated with lead 77: 

a) for normal lines: 1 if the line is going from the 
lower to the upper, —1 if it is going from the upper 
to the lower propagator, and otherwise; 

b) for anomalous lines: 1 for incoming lines within 
the upper and outgoing lines within the lower prop- 
agator, — 1 for for outgoing lines within the upper 
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and incoming lines within the lower propagator, 
and otherwise. 



The diagrammatic rules are formulated generally 
enough to account for any choice of chemical potentials of 
the leads, i.e., we allow for any bias voltages between any 
pair of leads. A variety of interesting phenomena, how- 
ever, shows up already when all superconducting leads 
are kept at the same chemical potential (set to per 
definition), and a noncquilibrium situation is generated 
only by applying voltages between the normal leads and 
the superconductors. It is this limit that we are going 
to analyze in the rest of the paper. In this case, the 
diagrammatic language simplifies further. As already in- 
dicated above, we simply can drop all information associ- 
ated with the Cooper-pair numbers in our diagrammatic 
rules. 



C. Green's functions 

If all superconducting leads are kept at the same chem- 
ical potential, which we put to per definition, we can 
use the simplified diagram, where the information about 
the Cooper pairs is ignored. Such a procedure is identical 
to having dropped the Cooper-pair states from the very 
beginning in the Hamiltonian, i.e., using the Bogoliubov 
transform without employing the operators . For this 
case, the charge current in lead r\ has been related to 
the local Green's functions of the quantum do t 14 i 24 ' 29 by 
using the approach of Ref. HH Here, we report the for- 
mula of Ref. which is useful for what comes in the 
following. For the sake of keeping the notation compact, 
we use the Nambu representation for the dot operators: 



given by 



(d h d\) 



The current flowing out of lead r\ is written as the sum 
of two contributions, J n = J\ n + J2^, with 



J. 



1,, 



1 1 



J211 — — 



(2G R H/,H + G<H)]} 
e ( duj 
h 2n 



: r r ,D I? (w)Re \ Tr 



(16) 
17) 



V" 2 - 
IAJ 



|A„|-M) : 



1 and 



for the superconducting leads, and D v (oj) 
D 7) (o;) = if t] describes a normal lead. 

The current J\ v involves only excitations energies u 
above the gap. This is the only contribution in a normal 
lead, where it reduces to the result presented in Ref. I33L 
For a superconducting lead, Ji^ has a contribution due to 
the normal elements of the dot Green's function, which 
describes quasiparticle transport and is independent of 
the superconducting phase difference, and a contribution 
due to the anomalous components of the Green's func- 
tions, which is in general phase dependent. 

On the other hand, Ji^ involves only excitations ener- 
gies lu below the gap and it describes both Josephson as 
well as Andreev tunneling. 

The above current formula becomes particularly use- 
ful in the limit of a large superconducting gap (|A^| — > 
00), where quasi-particle excitations are inaccessible, Ji^ 
dominates the transport. In this case, the current in the 
superconducting lead r\ reads 

J v = —TnKd^)] sin(*-$,), (18) 

where (did*) — \(d\dt)\ exp(i'I') is the dot pair ampli- 
tude. Equation (|18|) has a very simple meaning: it de- 
scribes the Josephson current between the lead with su- 
perconducting phase <i>^ and the dot with a phase ^. 
All the complicated physical effects due to the interplay 
of Coulomb interaction, coupling to all (normal and su- 
perconducting) leads and non-equilibrium due to a finite 
bias voltage between normal and superconducting leads, 
are hidden in the dot pair amplitude. 



III. RESULTS 

In the remaining part of the paper, we illustrate our 
formalism by considering two examples. 



where A v = y ^* ^QJ' anc ^ fvi 1 ^) = [1 + exp(w — 

MjjV^bT)] -1 is the Fermi function, with n v being the 
(electro-) chemical potential of lead 77 (= for the super- 
conductors), T the temperature and fee the Boltzmann 
constant. The local dot Green's functions G R (w) and 
G < (uj) are matrices in Nambu space, whose components 
(G < (w)) m and (G R (w)) are defined as the Fourier 

transforms of i(<j>l(<S)(t> m (t)) and -i6(t)({<t> m (t), <j>i (0)}), 
respectively. 

The two weighting functions D^oj) and D v (uj) are 



A. Josephson coupling due to cotunneling 

First, we analyze the equilibrium Josephson current 
through a superconductor-dot-superconductor system in 
the limit of weak tunnel coupling. The lowest-order 
mechanism that establishes a Josephson coupling be- 
tween the superconductors is cotunneling, i.e. the 
Josephson current starts in second order in the tunnel- 
coupling strengths T n . We consider a symmetric setup 
with both tunnel-coupling strengths equal to Ts and 
Al = = I A| exp(i<i>/2). The two superconductors 
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are kept at the same chemical potential /Ug = 0. We 
determine the Josephson current Jj OS = Jl = — Jr to 
second order in Ts, 



Jos 



i E [ 



W xx" L(2)p x'" (0) 



TrrXx"L(l)pX"(l) 

xx' x' 



(19) 

For the limit |A| ^> fee? 1 considered here, there is no 
microscopic mechanism in our model to make the dot 
degrees of freedom relax to equilibrium. This situation 
occurs because the quasiparticles excitation in the su- 
perconducting leads are not accessible. In reality, the 
degrees of freedom of the dot will be coupled to some 
thermal bath with temperature T and the dot will reach 
an equilibrium distribution. Hence, we assume that in 



zeroth order only the diagonal probabilities are non van- 
ishing and they are given by the Boltzmann factors 



P (0) = 

x 



cxp[-P x /(fc B T)] 



(20) 



with E = 0, P T = E i = e, E D = 2e + U and Z = 
J2 X ex P[ — E x /(k&T)\. Notice that in the model studied 
later in Section UlI B\ the presence of a normal lead tunnel 
coupled to the dot provides a mechanism for the dot to 
reach equilibrium. First, we focus on the regime that 
both e and e + U lie inside the superconducting gap. The 
only non-vanishing first-order correction to the reduced 
density matrix concerns the off-diagonal element Pq = 
(Pq 3 )* and it reads 



P 



0(1) 



2e + U 
2 

27+ U' 



IF, 



00(1) p(0) 



DO 







W 



0D(l)p(0) 



DD 



D 



r S COS - {A \rr-r )P 



>(0) 



.4 



■ U 



P 



(0) 



.4 



■ U 



p(0) 



(21) 



where the function A(z) is given by 

1 f°° 1 1 
A(z) = - / dx . (22) 

Equation (|21[) describes how a finite pair amplitude in 
the dot can be established in first-order in Ts- In fact, 

I 



(Pp)* = Pq = (d^df) is equal to the pair amplitude in 
the dot. 

Evaluating the second-order current diagrams 
l^xx' L *- 2 ' ) (an example is shown in Appendix [3J and 
using Eq. (TIT)]) , we obtain the following lengthy but 
complete result for the Josephson current 



JOS 



2e , 

— Tl sin$ 



:F - 



e + U 



——A 2 I — 
2e+U VIA' 



,(o) 



1 

Taj 



F - 



2e + U 

e + U 



-A 2 - 



U 



P 



(0) 



D 



4P 



— e e 
lAI'lA 



P 



p(0) 



(23) 



where the functions F(z) and B{z,z') are defined as In the limit |A| — > oo, all second-order current rates 

vanish and the Josephson current is given by Pq ^ mul- 
^ poo ^ tiplied by the corresponding first-order current rates, 

F(z) — — dx - - which yields 



Vx^~T 



1 1 1 1 ,„„ X 

rfy r^— — ; — ( 24a ) 



J jos =|rIsin$--^-(pr-pW). (25) 

i /-oo ^ The result Eq. (j2"3"|) for the second-order equilibrium 

B(z,z') — — / dx = Josephson current is valid when both the level e and e + U 

7T Ji Vs; — 1 are inside the gap, therefore the limit of large interaction 

1 }_ * (24b) U ^ co cannot be obtained directly from Eq. (j2"3")l . How- 

^/y 2 — l x + z x + y x + z' ever, in the limit of large interaction the double occupa- 
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FIG. 2: (Color online) Density plot of the Josephson current 
as a function of the level position e and of the interaction 
strength U. The region where the system behaves as a ir- 
junction is indicated by the symbol n in the plot. The dashed 
lines delimit the region where at zero temperature the dot is 
singly occupied. The other parameters used in the simulation 
are: fc B T/]A| = 0.05, r s /|A| = 0.01, $ = tt/2. 



tion of the dot is forbidden, and the Josephson current 
can be obtained by dropping all diagrams involving the 
doubly occupied state |D): 



J'ios — 



IT, 



00L(2) D (0) 



P n W + 2Wr L(2) P( 0) 



— sm$ 



F 



P, 



(o) 



F 



e 



Pi?2 3 



Equation (|26|) agrees with the results of Glazman and 
MatveevJi 

The Josephson current, Eq. (|2"5)l . is plotted in Fig. [2] 
as a function of gate voltage and interaction strength. 
We find, in agreement with Ref. y, the formation of a 
7r-state for gate voltages such that — U < e < 0, with the 
transitions being smeared out by temperature. 



B. Andreev-level spectroscopy 

We now turn our attention to the setup shown in Fig. [3] 
As compared to the geometry considered so far, there is 
a third, normal (N), lead with tunnel-coupling strength 
Tn, in addition to the two superconducting ones (L, R). 
Again, we assume the same tunnel coupling Tg and chem- 
ical potential /is = for both superconducting leads, and 
A L = = |A|exp(i$/2). The third lead allows for 
driving the quantum dot out of equilibrium by applying 
a voltage between normal and superconducting leads, ex- 
pressed by a non vanishing chemical potential /in of the 
normal lead. The quantities of interest are the Josephson 
current Jj OS = (Jl — </r)/2 and the Andreev current in 
the normal lead J an d = —{Jr + Jh)- 




R 



FIG. 3: Schematic setup of a quantum dot tunnel coupled 
to one normal and two superconducting leads. The dot can 
be driven out of equilibrium by a bias voltage applied to the 
normal lead. 



In Ref. [29| we studied this setup in the limit of weak 
tunnel couplings. We found that by applying a bias volt- 
age between normal and superconducting leads one can 
induce an out-of-equilibrium proximity effect in the quan- 
tum dot, which, in turn, supports a Josephson coupling 
carried by first-order tunnel processes instead of second 
order (cotunneling). We described the non-equilibrium 
Josephson current as well as transitions from to tt- 
states perturbatively to first order in Tg. This limited 
the applicability to a small range of gate voltages and 
temperatures larger than the tunnel-coupling strengths. 
The proximity effect was of purely non-equilibrium ori- 
gin since the influence of the superconducting leads on 
the quantum-dot spectrum, typically associated with the 
picture of Andreev bound states, could not be resolved. 

In the present paper, we want to go beyond the limit 
considered in Ref. [29| for two reasons. First, we aim at 
covering the full range of gate and bias voltages, thus, in- 
cluding both equilibrium and non-equilibrium proximity 
effect. Second, we are interested in mapping out the spec- 
trum of Andreev bound states of an interacting quantum 
dot. To pursue both of these aims, we need to go be- 
yond first-order transport in Tg . As usual for interacting 
systems, the full problem for arbitrary values of |A| and 
Fg can only be solved approximatively. In the limit of 
a large superconducting gap |A| — * oo, however, we are 
able to derive an exact result by resummation of the con- 
tributions of all orders in Tg . This is possible because for 
| A | — > oo only a small subset of all diagrams contributes 
to the generalized rates: the only superconducting lines 
that remain are anomalous ones that connect vertices 
within one (the upper or the lower) propagator, with no 
other vertex appearing in between. This simplification is 
related to the inaccessibility of quasiparticle excitations 
in the superconducting leads and the fact that a Cooper 
pair should tunnel in a time interval oc fi,/|A|, which be- 
comes infinitesimal for |A| — > oo. A rigorous proof is 
given in Appendix [Bl 
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In the limit |A| — > oo we can evaluate the current in 
the superconducting leads L, R by means of Eq. fig]) , i.e., 
we only need the pair amplitude of the quantum dot. It 
is useful to introduce a dot isospin defined as^ 



pD 



pO 
D . 



pD 

: 



pO 



Pn 



Po 



(27) 



Finite x- and y-components of the isospin indicate co- 
herent superpositions of the dot being empty or doubly- 
occupied. 

We rewrite the master equation for the dot reduced 
density matrix in the form of a Bloch equation for the 
isospin, taking into account all rates up to first order in 
Tn [the order in is indicated by the superscript (i) 
with i = 0, 1]. The Bloch equation for the isospin reads: 







dl 

It 



A-R I + Ix B 



(28) 



tive magnetic field acting on the isospin has a zeroth- 
order component B^ - 1 and a first-order component B^, 



(o) 



B 



B« = 



2r s cos$/2 \ 

(30) 
-(2e+U) J 

/-I m {< D(1) -< D(1) }\ 

. (31) 

v Im{< D «} ) 



The explicit expressions for the generation vector and 
the relaxation tensor can be written in a compact way, if 
we define the Andreev bound-state energies. These are 
given by the poles of the retarded Green's function of the 
dot for vanishing coupling to the normal lead, 



where the first, second, and third term describe genera- 
tion, relaxation, and rotation of the isospin, respectively. 
The explicit expressions of the needed generalized rates 
are reported in Appendix [D] In order to decouple the 
equations for the isospin from those for the diagonal prob- 
abilities we made use of the relations: W, 



D0(1) 



(1) 
Dct 



< D D(1) = o. 



The relaxation tensor and the generation vector start in 
first order in Tn ■ The generation vector reads 



A (D = 



Re{</ (1) } 



-Im I W, 



IT, 



(i) 



W, 



(i) 



OfT 



\ 



/ 



(29) 



The only non vanishing elements of the relaxation tensor 



are: R^l = 
and Rxz = 



B$ = -Re 



{< D(1) },4V 



= w. 



(1) 



■W. 



(1) 



p(l) 



ReiW, 



,0D(1) 



on 



DD(1)\ 

do y 



crO ~ vy itD ' 

The effec- 
I 



E A . 



l' , 7 




-) +rlcos2-, (32) 



where 7 and 7' can take the values ±1. There are four 
resonances which lie pairwise around zero energy. We get 
for the generation vector: 



Ai 1 



v 



Ai 1 



r s r N $ , ( „ , 

cos- 2^ 7/n(^A,7', 7 ) 



4e A 



7,7'=± 



= 




1 + 7 



7,7' =± 



(33a) 
(33b) 

/n(-Ea, 7 ', 7 ) ~ S 33c ) 



with 6a = y (e + + FgCOS 2 ^, where the square- 
root dependence clearly indicates that the result is non- 
perturbative in Tg. The non- vanishing elements of the 
relaxation tensor are: 



(i) _ 



±L ZZ 



a& = T E (1 

7,7'=± 



r N 



2 



E (1+7 



7,7'=± 



e + U/2 



■ cos „ 
2e A 2 



C-A 



e + U/2 



PA 



- 77n(-Ea,7',7) 



-77n(^a, 7 ', 7 ) 



■w E 777n(-Ea,7', 



'a) 



(34a) 
(34b) 
(34c) 



7,7'=± 



By means of Eq. (fT8| , the current in the su- component of the isospin provide the Josephson and An- 
perconducting leads can be written as Jr,l = 
(l y cos sin where the upper (lower) sign 

refers to the right (left) lead. Hence, the x- and y- 
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dreev currents, respectively, 

2e $ 

4e $ 

7 an d = -— r S / y COS— , (36) 

whereas the z-component is related to the charge in the 
quantum dot, 

Q = -e(l + 2/ z ). (37) 

We solve for the stationary solution for the isospin. 

Expanding the Eq. (|28|) to zeroth order in Tn yields = 
1(0) x B (0) ) andj thuSj j(o) || B (o)_ Tq determine the 

proportionality constant, we multiply B( ' from the left 
to Eq. ((28)) expanded to first order in Tn, and obtain the 
zeroth-order result 

1(0) - ( A(1) - B(0) ^ B(0) (38) 

which yields the Josephson current and the quantum- 
dot charge. The Andreev current, on the other hand, is 
proportional to the y-component of the isospin and starts 
in first order in T^. The first-order contribution to the 
y-component of the isospin can be derived by multiplying 
either x or z from the left to Eq. (f28|) expanded to first 
order 

tf> = -Lz.(A«-R«.I<°>) (39) 
= Lx-JaO-RW.^) . (40) 



The formation of a finite pair amplitude of the dot is 
favored if the empty and doubly-occupied dot states are 
degenerate, 2e + U = 0. In this case, however, the dot is 
preferably singly occupied in equilibrium, i.e., the prox- 
imity effect is strongly suppressed by Coulomb charging. 
For finite values of the superconducting gap |A|, a small 
Josephson current through the dot can be established by 
cotunneling processes. In the limit of infinite |A|, how- 
ever, this is not possible, and the proximity effect and, 
thus, the Josephson current is exponentially suppressed. 
In fact, we find for this regime = 0, i.e., no isospin 
is generated. 

The are two routes towards the generation of a finite 
dot pair amplitude. One is to change the gate voltage 
such that empty or double occupation of the dot becomes 
available. Then, the tunnel coupling to the supercon- 
ductors give rise to an equilibrium proximity effect that, 
however, starts in higher order in the tunnel coupling 
strength. To achieve a finite pair amplitude at lowest 
order already, one has to apply a finite bias voltage at 
the normal lead. This induces a non-equilibrium prox- 
imity effect that supports a first-order Josephson current 
through the dot. 

1. Equilibrium 



First, we consider the equilibrium situation (/xn = 0). 
In this case, the relation x • RW • B^/z • R^ • = 

Ax^ / Az ensures that no current flows in the normal 
lead. The exact result for the equilibrium Josephson cur- 
rent in zeroth order in Tn reads 



J jos = -Igsin$ 



+ U/2 



7,7' = ± 



[/(#A, 



7,7' — ± 7,7'— ± 



(41) 



where f(to) is the Fermi function with zero chemical po- 
tential. Notice that the only role played by the normal 
lead is to provide a mechanism for the electrons in the 
dot to reach equilibrium. Expanding Eq. (|4"T|) to second 
order in Ts we recover the result of Eq. (l2~5l) . 

In equilibrium, the pair amplitude of the quantum dot 
(did-\) for the symmetric setup and symmetric gauge is 
real, i.e. = ^ since I y = 0. In panel (a) of Fig.Q] 

we plot the pair amplitude as a function of the level posi- 
tion e for zero temperature. In particular we note that a 
quantum phase transition occurs if {7/2 > Ts| cos($/2)| 
at e = e± = -U/2± y/{U/2) 2 - f§ cos 2 ($/2). The values 
e± where the transition takes place depend on the tunnel 
coupling Tg and on the superconducting phase difference. 
The pair amplitude in first order Ts, see Eq. (|21l) . which 



gives rise to the second-order Josephson current, exhibits 
the phase transition at different values of the level posi- 
tion. The behavior of the pair amplitude can be under- 
stood by considering the Andreev-bound-state configu- 
ration in each region (see inset of panel (a) of Fig. [4|. 
Panel (b) of Fig. 2] shows the charge of the dot as a func- 
tion of the level position. In this case, the first-order Tg 
correction vanishes. The full result for the charge shows 
that due to proximity effect the charge on the dot is not 
always quantized^ 
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FIG. 4: Dot pair amplitude (dyd-\) (a) and charge (b) as a 
function of the level position e at zero temperature. Note that 
the pair amplitude in equilibrium for the symmetric setup is 
real. The dashed lines refer to the pair amplitude and charge 
calculated up to first-order in Fs; for the pair amplitude the 
Oth-order contribution vanishes, while for the charge the first- 
order one does. In the inset of panel a) we show a schematic 
picture of the Andreev bound state energies £a,v,7 m the 
three different regions, where the values 7' 7 are indicated 
next to the level. The sign of the contributions of the bound 
states _Ea,7',7 to (didf) is determined by the index 7. We 
have depicted the levels with 7 = + by a solid line and those 
with 7 = — by a dashed line. The other parameters used in 
the simulation are: Ts/U = 0.5, and $ = n/2. 



2. N on- equilibrium 

Next we turn our attention to the non-equilibrium sit- 
uation (/in 7^ 0). Applying a bias voltage to the normal 
lead produces a finite current in N, which is sustained 
by Andreev-reflection processes. We do not give here the 
explicit analytical expressions for the zeroth-order Jj os 
and the first-order J an d since they are rather lengthy. In- 
stead, in Figs. O [7| we plot the Josephson current, the 
Andreev current, and the dot charge as a function of the 
level position e and of the chemical potential of the nor- 
mal lead /in, for different values of the tunnel-coupling 
with the superconductor. 

In Fig. [Ha) one can see how the Josephson current 
can be controlled by the chemical potential of the nor- 
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FIG. 5: (Color online) Density plot of the Josephson cur- 
rent (a), the Andreev current (b), and the charge on the 
quantum dot (c) as a function of the level position e and 
of the chemical potential /xn- In panel (a), the region where 
the system behaves as a n -junction is indicated by the sym- 
bol 7T. The dashed lines map the Andreev bound states: 
/iN = Ej^ >rf i „. The other parameters used in the simula- 
tion are: k B T/U = 0.01, T N /U = 0.005, T s /U = 0.1, and 
$ = vr/2. 
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mal lead. There is, first, a broad region set by the charg- 
ing energy in which the Josephson current is suppressed. 
Second, a 7r-transition can be driven both by the trans- 
port voltage and by the gate voltage controlling the level 
position. For fixed /l«n the transition occurs at e = — J7/2, 
i.e. when the energy of the empty and double occupied 
dot are degenerate. We remark that this transition is 
slightly shifted when higher-order corrections to the ef- 
fective field are included. In fact, near the transition 
— 2e + U becomes small and hence needs to 
be taken into account. This has been done in Ref. 
in the weak-proximity limit. Panel (b) of Fig. [5] shows 
the Andreev current. It is largest at e = — U/2 outside 
the region where charging energy suppresses transport. 
Panel (c) of Fig. \5\ shows the dot charge. We find a pro- 
nounced feature around e = — U/2 that is associated with 
generating a y-component of the isospin by rotation out 
of the z-direction. 

In Figs. [6] and El the coupling to the superconduct- 
ing lead is stronger, comparable to the Coulomb in- 
teraction strength, and the term proportional to Tg in 
y/{e + U /2) 2 + Tg cos 2 (<E>/2) becomes more important, 
leading to a more pronounced splitting of the Andreev 
bound-state energies. We stress here that the current in 
the normal lead as a function of both gate and trans- 
port voltage maps the energies of the Andreev bound 
states in the dot. Therefore, measuring the current in 
the normal lead allows to perform a an Andreev-bound- 
state spectroscopy and, hence, to gather information on 
the superconducting correlations induced in the dot. 



IV. CONCLUSION 

We have presented a real-time diagrammatic trans- 
port theory for systems composed of interacting quan- 
tum dots coupled both to normal and superconduct- 
ing leads. First, we have applied this theory to study 
the Josephson current through a quantum-dot tunnel 
coupled to two superconductors in second order in the 
tunnel-coupling strengths. In particular, we have studied 
how a 7r-phase develops with increasing on-site Coulomb 
repulsion. Next, we have considered a quantum dot cou- 
pled to one normal and two superconducting leads, in 
the limit of large superconducting gap. In this regime, 
all orders in the tunnel- coupling strengths with the su- 
perconductors can be summed. This enabled us to in- 
vestigate the strong-proximity regime. In particular, we 
analyze the Josephson current and identify the parameter 
regions where the system behaves as 7r-junction; the ir- 
transition can be triggered both by the dot level position 
and the bias voltage. We find also that a spectroscopy of 
the Andreev bound states of the system can be realized 
by measuring the Josephson current between tho two su- 
perconductors, the Andreev current in the normal lead 
or the charge of the dot as a function of both the dot 
level-position and the bias voltage. 
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FIG. 6: (Color online) Density plot of the Josephson cur- 
rent (a), the Andreev current (b), and the charge on the 
quantum dot (c) as a function of the level position e and 
of the chemical potential ^n- In panel (a), the region where 
the system behaves as a n -junction is indicated by the sym- 
bol 7r. The dashed lines map the Andreev bound states: 
/in = -Ea.7',7- The other parameters used in the simula- 
tion are: ksT/U = 0.01, T N /U = 0.005, T s /U = 0.5, and 
$ = vr/2. 
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FIG. 8: Second-order diagrams contributing to the current 
rate< 0L(2) . 
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APPENDIX A: SECOND-ORDER, FINITE |A| 



FIG. 7: (Color online) Density plot of the Josephson cur- 
rent (a), the Andreev current (b), and the charge on the 
quantum dot (c) as a function of the level position e and 
of the chemical potential /in- In panel (a), the region where 
the system behaves as a n -junction is indicated by the sym- 
bol 7T. The dashed lines map the Andreev bound states: 
/in = -Ea.7',7- The other parameters used in the simulation 
are: k B T/U = 0.01, T N /U = 0.005, T s /U = 1, and $ = vr/2. 



In this Appendix we show, as an example, the cal- 
culation of the second-order current rate Wqq L ' . The 
second-order diagrams contributing to this rate are 
shown in Fig. [8l the signs have been assigned making 
use of Rule 6 in Section III B 31 Using the diagrammatic 
rules of Section III B 31 we get 
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| A | 2tt ^/uj 2 - [Ap 7|A| 2 ^ 



do/ 
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■ w + ecJ + u/w+e 



APPENDIX B: 



DERIVATION OF THE RULES 

FOR |A| -> oo 



Here, we give a rigorous proof of the rules which in the 
|A| — > oo limit allow us to greatly reduce the number of 
diagrams to be considered. 

Rule (i): No vertex should be considered between the 
two vertices of a line with a superconducting lead. 
Le us consider a diagram where a vertex v' exists between 
the two vertices of a superconducting line with energy oj. 
Let the vertex v' be associated with a line with energy 
lo' . According to the diagrammatic rules 2 and 3, this 
diagram contains the factor— ~±t~ ±u±u* D(lu)\A/lu\ 
which upon integration over to vanishes as 1/|A|. On 
the other hand, if no vertex is inserted between the two 
vertices of the superconducting line, the diagram con- 
tains the factor— ± ^ D(uj)\A/ui\ which remains finite 
upon integration over to. 

Rule (ii): No line with a superconductor joining the 
upper and lower propagator should be considered. 
Let us consider the diagram where a superconducting 
line is running from the upper to the lower propagator 
and the vertex on the upper propagator is on the left of 
the one on the lower propagator. In virtue of rule (i), 
the diagram with the two vertices swapped, i.e. with 
the vertex on the upper propagator being on the right of 
the one on the lower propagator, also exists. These two 
diagrams cancels each other for |A| — > oo. 

Rule (iii): No normal line with a superconductor should 
be considered. 

Let us consider a part of a diagram with a state |xu) 
running on the upper part of the Keldysh contour and 
with a state |xi) on the lower part. In virtue of the two 
previous rules, there are only four possible ways of insert- 
ing a normal line with a superconducting lead, which are 
schematically depicted in Fig. [9l In the large-gap limit, 
the diagrams corresponding to the possible insertion of a 
normal line, have the same absolute value. But the dia- 
grams arising from the insertion in the lower propagator 
(shown in the second line of Fig. [9]) have an opposite sign 
with respect to the ones in the upper propagator (first 
line of Fig. [5]). Finally, it easy to prove that for any |x u ) 
and there are, for our single-level model, only four 



possible insertion: two in the upper propagator and two 
in the lower propagator. Hence, the sum of all these dia- 
grams vanishes. To clarify this point, let us consider the 
exemplary case that |x u ) = |0) and = | T) : then the 
possible insertions are: (a) with intermediate state | f) 





X, 



X, 



(a) 



X, 



X. 



(b) 





FIG. 9: Possible insertion of a normal line in the upper or 
in the lower propagator. Note that these insertions are a 
part of a larger irreducible diagram and hence there are other 
lines (not shown) running from one side to the other of the 
insertion. 



or | I) for the upper propagator; (c) with intermediate 
state |0) and (d) with intermediate state |D) for the lower 
propagator. 



APPENDIX C: CALCULATION OF A 
GENERALIZED RATE IN THE |A| -> oo LIMIT 



Here, we show in one example how all contributions in 
Ts can be summed up. We consider the off-diagonal rate 
in first order in Tn and we add all contributions 
in Tg. In particular, only diagrams with an odd number 
of anomalous lines on the upper propagator contribute 
to this rate. The first two diagrams are shown in Fig.fTOl 
The contribution with 2n + 1 anomalous line reads 
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Summing up all terms we get 



< 0(1) = 2zr N r s cos| J 



r N r s $ 

= l cos — 

e A 2 



dui 

— f N (uj) — 

27T (uj 



1 



iO+){w -e + U - i0+) - (T s cos($/2))" 
1 1 



uj + U/2 -e A -iO+ ui + U/2 + e A - i0+ 




q t D a 2 Q q 3 
W W D 



ill 



^2 ^3 




FIG. 10: First two contributions to W™ (1) . 



with e A = vMe + C^/2) 2 + rg cos 2 ($/2). 



where the Andreev bound-state energies read 



-Ea,- 



7 




2 ' ' v V ' 2 

Some of the off-diagonal rates start in zeroth-order Tn- 



T2 COS 2 -. 

s 2 



In particular, we have 



,D0(0) 



(< D D(0) ) 



0D(0) 
DD 



= <o D(0) 
= (<o D(0) )* 



I^pj"' = irscos$/2. We also need the real part of 
the first-order corrections to these rates. Notice that the 
following relations hold Re{V ™ (1) } = Rc{tT°° (1) } = 

Re{< D «} = Rc{<r} and Re{<?«} = 

Rejl^} = Re{< D «} = Ref^}- The 
first-order corrections read 



APPENDIX D: GENERALIZED RATES TO ALL 
ORDERS IN T s 

In this Appendix we give the expression for the gener- 
alized rates, which are necessary to compute the Joseph- 
son current in zeroth-order Tn and the Andreev current 
in first order. The diagonal rates start in first-order Tn 
and they are given by 

= ^£(i + 7^)[1-/n(£a,-, 7 )] 

w% = ^2(i + 7 i±p)[i-^A, + , 7 )] 

O = ( 1 + 7 1 ^^)/N(i?A, + , 7 ), 



Re{< °«} = -r s cosf^]T 7 / N (£ A ,_, 7 ) 

7=± 

Re{< D D «} = -rscos|^^ 7 / N (^ A , +;7 ). 

7=± 



There are also some off-diagonal rates which start in first 
order in Fn: 
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Re{</«} = -Re{< D «} 



£sTn 
4e A 



cos it 2^ 7/n(-B A)7 ' i7 ) 



Re{< D «} = Re{<«} 



r N 



E 

7l7 '=± 



2 

1-7 



7.7'=± 

e + U/2 



7'/n(-Ea, 7 ', 7 ) - 
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